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Abstract 

We study analytically the highly damped quasinormal modes of Kerr black holes 
in the small angular momentum limit. To check the previous analytic calculations 
in the literature, which use a combination of radial and tortoise coordinates, we 
reproduce all the results using the radial coordinate only. According to the earlier 
calculations, the real part of the highly damped quasinormal mode frequency of Kerr 
black holes approaches zero in the limit where the angular momentum goes to zero. 
This result is not consistent with the Schwarzschild limit where the real part of the 
highly damped quasinormal mode frequency is equal to c 3 ln(3)/(87rGM). In this 
paper, our calculations suggest that the highly damped quasinormal modes of Kerr 
black holes in the zero angular momentum limit make a continuous transition from 
the Kerr value to the Schwarzschild value. We explore the nature of this transition 
using a combination of analytical and numerical techniques. Finally, we calculate 
the highly damped quasinormal modes of the extremal case in which the topology of 
Stokes/anti-Stokes lines takes a different form. 



1 Introduction 



Black hole quasinormal modes (QNMs) are the natural vibrational modes of perturba- 
tions in the spacetime exterior to an event horizon. QNM frequency spectrum is composed 
of an infinite number of discrete and complex frequencies, u n = ujr + ic^j, with the "over- 
tone" number n — 1, 2, 3, The imaginary part of the frequency, ui, signals the presence 

of damping, a necessary consequence of boundary conditions that require energy to be 
carried away from the system. 

Black holes are used for exploring ideas in quantum gravity, and their QNM frequency 
spectrum is an obvious place to search for a quantum signature. One possible link between 
QNMs with high damping rates and the semi-classical level spacing in the black hole 
quantum area spectrum has been proposed originally by Hod[lJ in 1998 and modified 
recently by Maggiore[2]. Many people have since calculated the highly damped QNMs of 
different black hole models. 

The QNMs of Kerr black holes with high damping rates were first explored numerically 
by Berti et al. in [31 HJ E]- It was not until recently that Keshet and Hod[6j were able to 
analytically calculate the highly damped QNMs of Kerr black holes for the first time. This 
work followed with a more detailed paper by Keshet and NeitzkeJT]. Kao and Tomino[8] 
generalized the results in [6] to spacetime dimensions greater than four for scalar (spin-zero) 
perturbations. According to the results obtained in [31 HI El El E] , the real part of the 
highly damped QNM frequencies of Kerr black holes approaches zero in the Schwarzschild 
limit of a — > 0, where a is the black hole angular momentum per unit mass. This is in 
clear contradiction with the Schwarzschild result in which the real part of QNM frequencies 
approaches ln(3) / (8nM), where M is the black hole mass, as the damping rate approaches 
infinity. (Here, we use geometrized units where G — c — ks — 1.) In other words, 
it appears that the large damping limit (\uii\ — > oo) and the zero angular momentum 
limit (a — > 0) of Kerr QNMs do not commute. In [6], Keshet and Hod state that "the 
asymptotic QNMs are not continuous at a = 0" . In other words, at a = the real part of 
the QNM frequency makes a discontinuous transition from zero to the Schwarzschild value 
of ln(3)/(87rM). This issue was presented by Maggiore in [2] as an argument against Hod's 
proposal in [1] and consequently provides one of the motivations for this work. Note that 
the imaginary part of the asymptotic Kerr QNMs does not suffer from such discontinuity. 

Musiri and Siopsis[9] show analytically that when the range of QNM frequencies is 
bounded from above by 1/a (\au\ < 1), the Kerr QNMs do coincide with the Schwarzschild 
QNMs in the high damping limit. The calculation in [9] is valid for small values of the 
parameter a and includes the Schwarzschild case where a = 0. According to [9], in the 
large damping limit 

ln ( 3 ) 

where m is the azimuthal eigenvalue of the wave. 

In this paper we use a different analytic technique to reproduce the result obtained 
by Musiri and Siopsis in [9|. We also reproduce the result obtained by Keshet and Hod 
in [B] using a slightly different analytic technique in which the tortoise coordinate is not 
required. We show that the result obtained in [SI U\ is only valid when \au\ > 1. We also 
show that for any large but finite QNM frequency, the transition from the Kerr value to 
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the Schwarzschild value happens in a continuous way. We explore in detail this transition 
which happens at around \au\ ~ 1. Such transition becomes discontinuous only when 
\oj\ = oo. In this paper, we also calculate the highly damped QNMs in the interesting case 
of extremal Kerr black holes for the first time. 

In section 2, we set up the problem. We present the analytic calculations in section 3. 
In section 4, we explain the numerical calculations and provide the results. We conclude 
the paper in section 5. 



2 Wave Equation 



Following the notations used by Keshet and Hod in the radial part of Teukolsky's 
wave equation for a Kerr black hole with mass M and angular momentum J can be written 
in the general form 

+ P(r)tt(r) = , (2) 



where 



We define 



P(r) 



dr 2 

_ qo{r)oJ 2 + q\{r)u + q 2 {r) 
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(3) 
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(r 2 + a 2 ) 2 



qi 



2am(2Mr - q 2 ) + 2is r(A + q 2 ) - M(r 2 - a 2 ) 



(4) 
(5) 
(6) 

q 2 = m 2 a 2 - (A lm + s)A + M 2 - a 2 - q 2 - s(M - r) [2iam + s(M - r)] , (7) 

where a = J/M, I and m are angular and azimuthal harmonic indices, A[ m is the separation 
constant, and s — 0, —1/2, — 1, —2 is the spin-weight parameter for scalar, two-component 
neutrino, electromagnetic, and gravitational fields respectively. Note that Teukolsky's 
wave equation can be generalized to include black holes with electric charge q[Hj\, where 
the case of q ^ in the above equations is understood only for scalar fields. The roots of A 
determine the radii of outer and inner horizons. These roots are r± = M±y/M 2 — a 
The differential equation ^ has two linearly independent WKB solutions: 

< 



q- 



where 



P(r) 



(9) 



4(r — r_) 2 4(r — r + ) 2 

The last two terms on the right hand side of the above equation come from matching the 
WKB solutions to the exact solutions of the wave equation ^ in the limit where r — >■ r±. 

The boundary conditions that we apply to the wave equation ^ are outgoing-wave at 
infinity and ingoing-wave at the horizon. We will choose the phase of the square-root of 
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Q 2 such that ^fQ 2 ~ w as r oo. This means that the outgoing- wave solution at infinity 
is proportional to f\ while the ingoing- wave solution at r + is proportional to / 2 . Here, we 
assume that perturbations depend on time as e~ tu)t . This means that uj < 0. 

In the large damping limit (u>i — > — oo), the angular separation constant can be ap- 
proximated to be 

A lm = iA lC LUJ + (A + m 2 ) + O{\uj\- 1 ) (10) 
and in the Schwarzschild limit where a — > 0, this constant can be approximated to be 

A lm = 1(1 + 1) - s(s + 1) + 0(au) . (11) 

It is easy to show that in the limit a — > 0, and as long as \u\ <C 1/a, we get 

q w r 4 , 

gi ~ 2?s(r 3 - 3Mr 2 ), 

1 - s 2 

g 2 ~ -A Zm (r 2 -2Mr) + ^^(4r 2 + 4M 2 -8rM)-(l-s)(r 2 -2Mr) . (12) 

If, in addition to the above conditions, we assume that \u\ ^> 1, we can simplify go, Qi, and 
q 2 in the expression q u 2 + q±uj + q 2 further to find that 

go ~ r A , 
g 1 ~ —QisMr 2 , 
1 - s 2 

g 2 ps ^^(2M) 2 . (13) 

Note that we cannot ignore the term qiu because this term becomes comparable to gow 2 
and g 2 when r is of the order of magnitude of |w| -1 / 2 . 



3 Analytic Calculations 

In order to extract the WKB condition on the asymptotic QNM frequencies, first, we 
need to determine the zeros and poles of the function Q 2 and consequently the behavior of 
the Stokes and anti-Stokes lines in the complex r-plane. Stokes lines are the lines on which 

Re £ Q(r')dr' = , (14) 

and anti-Stokes lines are the lines on which 

Im j r Q{r')dr' = , (15) 

where t is a zero of the function Q 2 . The poles are located at r±. 
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Figure 1: A schematic illustration of Stokes (dashed) and anti-Stokes (solid) lines for finite a. The hollow 
circles represent the zeros of Q 2 , while the filled circles are the two poles. The thin line is the path taken 
along the anti-Stokes lines to derive the WKB condition on uj. 



3.1 The case of finite a: 

In this case the zeros of the function Q 2 are, by a factor of close to the zeros of 

the function go- After determining the topology of the Stokes/anti-Stokes lines, we follow 
the path shown in Fig. [T] starting on an unbounded anti-Stokes line at a. The solution on 
line a is known due to the boundary condition at infinity: 



ft 



(*a) 



(16) 



We apply the rules that are explained by Andersson and Howls in [TT] to move along the 
anti-Stokes lines. The steps we take are as follows: 



/r - in 
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(18) 
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where t 

712 = -721 = / 2 Qdr (22) 
Jt 

along a path to the left of the event horizon and 



ti 



7i2 = r Qdr (23) 
Jti 

along a path to the right of the event horizon. The boundary condition at infinity requires 
that the coefficient of f% in fy a be equal to zero. This leads to the WKB condition 

e 2ni2 = _ lm (24) 

Note that the above condition gives us the correct monodromy as a result of rotating 
clockwise around the outer horizon in the complex r-plane, which is 

= e~ iT V a , (25) 

where 

T = 712 + 721= j> Qdr. (26) 

clockwise 

In the large cu limit, we can write 

2«7i2 ~ 2iu / ^^—dr + 2i / ——rdr = uS Q + is5 s + iAxSa + rn5 m , (27) 
Jti A Jti 2^/q n l\ 



where 



4 = 2i f"A r , Sm ^ 2i p =^MLzm dr , 

Jt! A Jt x 



f* 2 -a , , p /■*! r(A + Q 2 - M(r 2 - a 2 ) , 

6 A = 2i dr, and 5 S = 2i V — ; — ^ -dr . (21 

Jt! 2J% Jti V^o A 



It is now easy to show that the WKB condition (24) results in 

oj = -my — i(<j) — n5) , (29) 
where n is a large integer and Co = 5 m /8o, 5 = 2n/5o, and <j> = (s5 s + Ai5a — 7r)/5o- Equation 



(29) is the same equation that Keshet and Hod found in [6]. 



3.2 The case of a — > 0: 

In the limit a — >■ 0, as long as >> 1 and \tua\ « 1, we can combine equations 
Q, and to see that the four zeros of the function Q 2 approach 



± / (3 ± 2V2)M S 
V c«j 
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Figure 2: A schematic illustration of Stokes (dashed) and anti-Stokes (solid) lines when a — > 0. The 
hollow circles represent the zeros of Q 2 , while the filled circles are the two poles. The thin line is the path 
taken to derive the WKB condition on u. 



Since s < for all types of perturbations and uj ~ — the four zeros above are located 
on the real axis. The structure of Stokes and anti-Stokes lines are shown in Fig. [2j In order 
to derive a WKB condition on the QNM frequency u we will follow the path indicated by 
the thin line shown in Fig. [2] along anti-Stokes lines. We begin our path on the anti-Stokes 
line labeled a which extends to infinity. According to the boundary condition, the solution 
on this line is 

*a = ■ (31) 

The steps that we take are very similar to the the steps explained in the previous subsec- 
tion. We choose our branch cuts in a way that the WKB solution fx is dominant on the 
Stokes lines extending to infinity and f'2 is dominant on the Stokes line ending at the outer 
horizon. After making a complete loop around the horizon, we return to the anti-Stokes 
line a. The final solution that we get is 

\]>_ _ e _i 744 ^ _|_ e -2i7i2 e -2i723 e -2i734 _|_ e -2i723g-2«734 _|_ g-2i7 3 4 _|_ e 2i 744^j 
_j_^ e -«744 A _|_ g2i7i2 _j_ e 2i7i2 e 2i723 _|_ ^712^1723^734^ 

^-2*712 g-2i723 g-21734 _|_ e -2i723g-2«734 _|_ e -2«7 3 4 _|_ e 2i 744"j ^32) 



The boundary condition at infinity requires that the coefficient of f2 in solution (32) be 
equal to zero. This leads us to the WKB condition 

g2i744 _ _ e -2i7i2 e -2i723g-2i734 _ e -2«723g-2i734 _ e -2i7 34 ^33^ 



where 



744= j Qdr = -2m Res Q = -AniuM (34) 

clockwise 



and in the large damping limit 




QisMu s 2 M 2 



1/2 



c 2 -=l-^) dr. (35) 
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Figure 3: A schematic illustration of Stokes (dashed) and anti-Stokes (solid) lines for the extremal case. 
The hollow circles represent the zeros of Q 2 , while the filled circle is the pole at the point where the inner 
and outer horizons coincide. The thin line is the path taken along the anti-Stokes lines to derive the WKB 
condition on w. 



Note that condition (33) gives us the correct monodromy as a result of rotating around 
the outer horizon in the clockwise direction, which is 



(36) 



The integral (35 ) can be solved analytically. We use the change of variable y = lot 2 / sM 
to get 



S f(3±*V2)i ( Qi i 

y y 2 



1/2 



4 J(^2y/2)i 



dy 



y 2 — 6iy — 1 — arccot [3y — Jy 2 — 6iy — 1 — 3iarcsinh 



-3i + y x 
2y/2 , 



+ 1 ln / y 

2 \3i — 5y + 3y/y 2 — Qiy — 1 



-| (3^:2^)4 



- (3^2^)4 



S 

^2 



Combining Eqs. (33) and (37), we find 

e 2i744 = _ 3 

for gravitational and scalar perturbations and 

e 2i744 _ _l 



(37) 



(3f 



(39) 



for electromagnetic perturbations. These WKB conditions are in perfect agreement with 
the Schwarzschild case. 
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Figure 4: Numerically generated anti-Stokes (left) and Stokes (right) lines for the extremal case. The 
small hollow circle in the middle is the pole at the inner horizon and the bigger hollow circles are the 
location of the zeros of Q 2 . 



3.3 The extremal case: 



In the case of extremal Kerr black hole spacetime where the inner and outer horizons 
coincide (a = M), the topology of Stokes/anti- Stokes lines takes a new form shown in Fig. 
[3j In order to determine the WKB condition, we follow the path shown in Fig. [3j which 
begins on the anti-Stokes line labeled a that extends to infinity and after moving from t 2 
to t\ along the anti-Stokes line to the right of the horizon the path ends on an anti-Stokes 
line which connects to the horizon r^. The steps that we take will be identical to the ones 
we took in subsection 3.1 from Eq. (16) to (19). The boundary condition at the horizon 
requires that the coefficient of f\ in Eq. ( 19 ) be equal to zero. This will lead us to the same 
WKB condition that we found in Eq. (24). This shows that the highly damped QNMs 
of Kerr black holes in the limit a — > M does indeed coincide with the highly damped 
QNMs of extremal Kerr black holes, even though the topology of Stokes/anti-Stokes lines 
in the extremal case is different. In other words, to find the asymptotic QNM frequency 
of extremal Kerr black hole all we need to do is to replace the angular momentum a with 
mass M in Eq. (29). This result is consistent with the high overtone QNMs of other types 



of extremal black holes [12 



4 Numerical Results 

The configuration of the anti-Stokes lines was found using the LinelntegralConvolutionPlot 
in Mathematica applied to the vector field (cos [at (r)], — sin[a(r)]) where a(r) = arg[Q(r)]. 
This vector field defines a direction for dr at every point in the complex plane such that 
aig[Q(r)dr] = 0. The anti-Stokes lines follow this vector field to and from the zeros. Sim- 




Iog[a,| 



Figure 5: The real part of the highly damped QNM frequency as a function of a (left) and the damping 
rate (right). Logarithmic scale is used along the horizontal axis to better illustrate the transition region 
from the Kerr value to the Schwarzschild value. The dots are the numerically generated data points and 
the solid line is a cubic spline displayed to suggest the possible behavior in the transition region. The 
hollow circle on the vertical axis indicates the location of the Schwarzschild value lor = ln(3) /8%M. Note 
that wr does not appear to be discontinuous at a = 0. We take M — 1/2, m = —1, ujj — — 10 6 (left), and 
a = 1CT 6 (right). 



ilarly, the Stokes lines follow the vector field (sin[a(r)], cos[a(r)]). Fig. [4] shows the vector 
fields for the anti-Stokes and Stokes lines for the extremal case illustrated schematically in 
Fig.|3j 

The numerical results in Fig. [5] use numerical integration to find the integrals 7y = 
Su Q{r)dr along a path in the complex plane homotopic to the anti-Stokes lines. We then 
use the WKB conditions (24) and p3j ) to calculate ujr in different regions. Thus the 
numerical techniques, as well as the analytical techniques, used in this paper are not valid 
when the anti-Stokes lines do not connect the zeros of Q{r). As the zeros move from the 
configuration seen in Fig. [I] to that seen in Fig. [2] the anti-Stokes lines no longer connect 
the zeros as can be seen by examining the vector field in Fig. [6j 

The results in Fig. 5 show that when a <C r^r, ujr —> ln(3)/(87rM) and for a ^> rjjr, 
the results agree with the numerical and analytical results in [31 HJ El El |8] . 



5 Conclusions 

In this paper, we used a different analytic technique to reproduce the result obtained by 
Musiri and Siopsis in In addition, in order to check the previous analytic calculations 
of the highly damped QNMs of Kerr black holes U\, we reproduced the results without 
the use of the tortoise coordinate. Our calculations suggest that, contrary to the previous 
results in the literature, the highly damped QNMs of Kerr black holes in the zero angular 
momentum limit make a continuous transition to the Schwarzschild value. See Fig. [5j 
A similar situation exists also in the highly damped QNMs of Reissner-Nordstrom black 
holes. In that case as well, it was initially thought that if we first consider the asymptotic 
QNMs of a Reissner-Nordstrom black hole and then let the charge go to zero, ujr does not 
reduce to ln(3)/(87rM), but rather reduces to ln(5)/(87rM)[13j. Such discrepancies in Kerr 
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Figure 6: Numerically generated anti-Stokes lines in the transition region from the Kerr topology to the 
Schwarzschild topology. The small hollow circle in the middle is the pole at the inner horizon and the 
bigger hollow circles are the location of the zeros of Q 2 . 

and Reissner-Nordstrom black holes were used by Maggiore[2j as an argument against the 
validity of Hod's conjecture [I]. However, similar to what we did here for Kerr black holes, 
it is shown in [T4] that the highly damped QNMs of Reissner-Nordstrom black holes in the 
zero charge limit do coincide with the known Schwarzschild results. 

If we try to derive the horizon area spectrum of a slowly rotating black hole assuming 
that the real part of the highly damped QNM frequency approaches zero when a — > 0, we 
will end up with zero spacing between the consecutive states of the spectrum according 
to Hod's interpretation of the highly damped QNMspQ. For that reason, Vagenas[16] 
and Medved [T7] use the imaginary part of the QNM frequency, which approaches the 
Schwarzschild value when a — > 0, in deriving the area spectrum of slowly rotating black 
holes following Maggiore's interpretation [2J. The results of this paper show that Hod's 
interpretation also gives a non-zero spacing between the consecutive states of the area 
spectrum. Following the steps of Medved in [T7] , we can replace u in the adiabatic invariant 
quantity 



where Q = a/{r 2 + + a 2 ), with ln(3)/(87rM) and we obtain an area spectrum for slowly 
rotating black holes of the form 




(40) 



A k + 0{4) «4Hn(3)Z| 



(41) 
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where k is a large integer. In the case of finite a neither Hod's interpretation nor Maggiore's 
interpretation of the highly damped QNMs give a unique value for the spacing in the area 
spectrum because both ojr and the spacing in uj are dependent on the parameter a. 

It is relevant to mention that Babb et al. in [TS] show that, when one includes quantum 
corrections to a Schwarzschild spacetime, ln(3) does not appear in the region of the QNM 
spectrum where the damping is larger than the inverse Planck/polymerization length scale, 
but it does appear in lower damping rates. They conclude that it is still plausible to 
think that the QNMs with damping rates that are large relative to the inverse horizon 
length M -1 ) but small compared to other inverse lengths will provide information 
about the black hole horizon. A similar conclusion can be made regarding the Kerr case. 
The QNMs with a damping rate much larger than M -1 but much smaller than a -1 may 
provide information about the mass of the black hole only. The information on the angular 
momentum of the black hole appears in the QNM spectrum of slowly rotating black holes 
when the damping rate becomes comparable to a" 1 . In other words, it seems that for Hod's 
conjecture to be true, we need to associate ln(3) with the mass of the black hole only. When 
a becomes comparable to the mass of the Kerr black hole, the two regions of the QNM 
spectrum overlap which makes it difficult to separate the information regarding the mass 
and angular momentum. Now, the question is, how can we relate the QNM frequencies in 
the intermediate damping region M _1 <C <C a -1 to the area spectrum of Kerr black 
holes if they only contain information on M. One way is to follow the intriguing suggestion 
of Makela et al. in [181 IIH] where the authors argue that in the case of multi-horizon black 
holes, it is the area of both inner and outer horizons which should be quantized in equal 
steps. For the Kerr case this area is 

Aotai = 4vr(r2 + a 2 ) + 4vr(r 2 + a 2 ) = WttM 2 , (42) 

which is the same as the area of a Schwarzschild black hole with mass M. We now can 
use the value of ur in the intermediate damping region of the Kerr QNM spectrum to find 
that 

AA total = 32vrMAM = 32vrM^^ = 4 ln(3)^ . (43) 

This is the same spacing that one finds in the Schwarzschild black hole area spectrum 
according to Hod's interpretation. Also, Maggiore's interpretation can only give a unique 
value for the spacing in the area spectrum if we use ui in the intermediate damping region, 
which gives 

AA total = 327rMh(\u I \ n - |w/| n _i) = Z2nM-^ = 8^ . (44) 

This result is the same as the result that one finds for the Schwarzschild case using Mag- 
giore's interpretation. 
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